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Abstract 

We describe how complex fluctuations of the local environment of an optically active quan- 
tum dot can leave rich fingerprints in its emission spectrum. A new feature, termed "Fluctua- 
tion Induced Luminescence" (FIL), is observed to arise from extremely rare fluctuation events 
that have a dramatic impact on the response of the system — so called "black swan" events. 
A quantum dissipative master equation formalism is developed to describe this effect phe- 
nomenologically. Experiments performed on single quantum dots subject to electrical noise 
show excellent agreement with our theory, producing the characteristic FIL sidebands. 
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Discrete solid-state quantum emitters such as III-V quantum dots^ and NV-centers' repre- 
sent highly versatile hardware for all-optical quantum technologies. In particular, self-assembled 
InGaAs Quantum Dots (QDs) are robust, ultra narrow linewidth sources of triggered single pho- 
tons and polarization-entangled photon pairs P The creation of remote entanglement via photon- 
interference'^ calls for discrete emitters that produce Fourier-transform limited single photon wave 
packets. This is very difficult to achieve in solid-state systems since the system typically under- 
goes frequency jitter due to fluctuations of the charge environment of the systemic and unavoid- 
able coupling to the lattice.'^ However, for low temperatures and weak optical driving, coupling 
to the lattice is inhibited^ and the spectral lineshape is expected to reflect slow fluctuations of 
the environment which manifest themselves as discrete spectral jumps and wandering.!^ Fourier 
transform-limited lines are typically not achieved in optical experiments, with measured linewidths 
somewhat above the theoretical limit.'^^^ While spectral fluctuations in self-assembled QDs have 
been investigated with non-resonant excitationj ^^*^^ ' the impact of a fluctuating environment in 
the case of true resonant excitation has received significantly less attention. Spectral fluctuations 
arising from environmental fluctuations are a common feature in condensed matter systems aris- 
ing also for colour centers in diamond,'^ colloidal nanocrystals^ and semiconductor nanowire 
quantum dots.f^ 

Here, we probe the impact of environmental fluctuations on the optical response of a resonantly 
driven quantum emitter embedded within an electrically tunable solid-state environment .1 ^^ * ^ ^ I The 
resonant excitation simplifies the picture since it avoids complications due to incoherently injected 
carriers, such as fluctuations arising from charge (de) trapping at defects in the wetting later. It 
also allows to consider only a single projection of spin by using circularly polarized excitation. As 
depicted schematically in[TJa), we consider a situation where a narrowband single frequency laser 
coherently excites the interband p-transition in a single dot and the emission from the s-shell is 
monitored.l^^lThe corresponding Hamiltonian reads: 

H = (A-^2i) |1) (1| +A|2) (2| +^2(|0) (2| + |2) (0|) (1) 
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where A = 0)2 — a)L is the detuning between the state 2 (interband p-transition) and the laser, 
and £"21 = 0)2 — 0)1 the energy gap between the p and s-interband transitions (we set /i = 1 for 
convenience). The relaxation between the p and s-shells — due to, e.g., phonon mediated pro- 
cessesf^S — occurs at the rate 721 and is included as an incoherent relaxation in the Lindblad form, 
T'^i)(2| where, for any operator c, is the Liouville superoperator defined as ^c(p) = 2cpc^ — 
c^cp — c^cp. The radiative decay of the s-state giving rise to luminescence is also described as 
where c = |0) (1|. The photoluminescence (PL) spectrum 5(0), A) is obtained by calcu- 
lating g(^)(t) = (cj'I'(O)cj(t)) from the master equation dtp = i[p,H] + ^^i)(2|(p) + ^J^ct(p) 
and taking the Fourier transform. Setting as the reference energy, the energy at which the s-state 
emits when A = (with maximum gain), the result reads: 



5(a),A) 



n (r^o + - A)2)(r2 + 4ft)2) + 8Q2(yiQr - 4(0) -A)(o) + 16Q4 



(2) 



where we introduced F = 710 + T^i. 
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Figure 1: (a) The level structure under study, (b) dressing of the p-state at resonance and (c) the 
MoUow doublet which arises from it. The laser is tuned closed to the excited p orbital, labelled 2. 
Incoherent relaxation 721 transfers the excitation to the s-state, labelled 1, which PL is detected. 
When the excitation rate is large enough, the laser dresses the states which are driven, namely, 
the p-state and the ground state. The latter, being common with the s-state, result in the observation 
of a MoUow doublet with splitting increasing linearly with Q.. 



This is the full quantum picture which describes saturation of the two-level system and dressing 
of the p-state by the laser when the pumping is large."^ In this case, this results in a spectral 
doublet rather than the MoUow triplet for a two-level system,'2Sl because relaxation is from the 
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naked s-state to the dressed ground state, which is common to both the bare s and driven p levels, as 
shown in[TJb). Provided that the gain is large enough, this configuration might have advantages as 
compared to resonance fluorescence, since excitation and detection are detuned and pump photons 
can be filtered from the emission. 

In the linear regime, Eq. (??) assumes the much simpler form: 

5(„,A) = ?5^xi^^|!^xi— J"^— j + 0(n)\ (3) 

no 7r(72i/2)2 + A2 ;r(7io/2)2 + (a)-A)2 ^ ^ ' 

with 0{Q!^) terms of at least fourth order in Q.. This shows that the coherently excited QD lumines- 
cence is given simply by the product of /) the emitted intensity, //) the population gain of the driven 
detuned oscillator, and ///) the Lorentzian PL lineshape of emission; an extremely simple and fun- 
damental result. The emitted intensity grows linearly with Q.^ as long as the system remains in the 
linear regime under its effective pumping Pgff. that is, as long as Gain(A) = Peff/7io <C 1. At reso- 
nance A = 0, where the driving is most efficient, the condition reads Gain(O) = (721710) <^ 1- 
We will restrict our discussion to the linear regime in the remainder of the manuscript, where one 
expects only emission from the s-shell directly to the ground state by direct radiative recombi- 
nation. Plotting the intensity of emission of the s-state photons (in colour, with lighter shades 
corresponding to higher intensities) detected at a given energy (on the x axis) for the various de- 
tunings of the lasers with the p-state (on the y axis), this produces a tilted line, since the emission 
linearly tracks the energy detuning from resonance. This feature is presented by the diagonally 
shifting transition shown in [2]^ a). The horizontal line superimposed at A = arises from the strong 
enhancement of the emission signal at resonance, termed "gain" in the discussion below. 

Measurements were made on many different single quantum dots and QD-Molecules (QDMs)'^^'^ 
to test the predictions of this theory. As discussed below and presented in[2fb), the experimental 
findings differ significantly from these simple expectations. The samples investigated were GaAs 
n-i-Schottky photodiodes containing a low density layer of InGaAs self-assembled QDs or ver- 
tically stacked QDMs. Such structures facilitate control of the static electric field in the vicinity 
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of the nanostructure via an applied voltage. The growth conditions used gave rise to QD nanos- 
tructures with a large radius (30-50nm) and, thus, a large oscillator strength as evidenced by mea- 
surements of short radiative lifetimes for self-assembled nanostructures (< 700ps, limited by the 
temporal resolution of our setup). Individual QDs or QDMs were excited optically using a tunable 
single frequency laser and their photoluminescence was recorded using a low temperature confo- 
cal microscope. Photoluminescence excitation (PLE) spectra from the same dot were obtained by 
recording the emission in the vicinity of s-shell transition with a CCD camera while scanning the 
laser energy through the excited orbital states of the system — fully analogous to the theoretical 
scheme introduced in the discussion related to[TJ A typical experimental result is presented in[2]^b) 
in a false colour representation for excitation in the range 9-13meV above the s-shell neutral exci- 
ton transition. An unexpected vertical line is clearly observed, labelled FIL in|2]^b), that is pinned 
to zero detuning. The energy gap between this newly observed FIL feature and the s-shell emission 
is identical to the detuning of the laser from the p- shell transition. In a configuration where the 
single mode laser was fixed and the p and s transition were shifted, for instance by application 
of an electric field to tune the levels through the DC Stark effect, we obtained identical results. 
This observation was also found to hold for QD-molecules pumping the antibonding branch of the 
neutral, coupled excitonl^SI and detecting the bonding branch. Here, the energy gap between the 
two transitions involved varies with electric but the FIL peak always appears at the energy 

of maximum gain. We will now show that this striking phenomenology is fully accounted for by 
fluctuations, albeit of a particular type. 

We proceed with the introduction of fluctuations in the theoretical model, an undertaking that is 
nontrivial at a microscopic leveL'^^lQur modelling remains rooted in the master equation approach 
and hinges on a reservoir that can be in any of Nr + I configurational states, labelled as /? = 
—Nr/2, . . . ,0, . . . ,Nr/2, each of which has the effect of shifting the levels of our system according 
to cOiR = C0i + R5, / = 1,2. Such a picture has been introduced by Budini^^ with a small reservoir 
to describe blinking. Here, we will assume that the reservoir is large enough to approximate a 
continuous fluctuating bath, introducing the need of a continuous distribution for the fluctuations. 
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Figure 2: (a) The expected luminescence from a resonantly driven QD in the linear regime: the 
diagonal transition shows the radiative recombination of the s shell as the frequency of the driving 
field is detuned from the excited state, while the horizontal line reflects the high increase of the 
excitation efficiency of the driven oscillator, (b) A new peak — FIL — appears in the experimental 
observation as a vertical line pinned by the laser, (c) The theory supplemented with fluctuations 
reproduces the observed phenomenology, (d) Scheme of the processes which give rise to FIL. The 
levels 1 and 2 have scalefree fluctuations which, when they bring the system into resonance, result 
in a strong absorption and subsequent emission. 



The driven QD Hamiltonian now depends on the states of the reservoir: 



Hr = (A -^21 +R5)\IR) {IR\ + {^ + R5) \2R) {2R\+^{\0R) {2R\ + \2R) (07?|) (4) 



where we use a symbol to denote that a fluctuating environment has been included in the picture. 
The states with / = 0, 1 or 2, correspond to the case where the QD is in the state |/) and the 
reservoir in its configuration R. Note that the latter is not a proper quantum number, but a label to 
index a macrostate for the environment. We assume that only the energy of the levels fluctuates 
while other parameters remain unchanged. In this case, the master equation of the system now 
reads: 

Tin y?! 
dtPR =i[pR,HR] + ^^\or){ir\{Pr) + Y^\iR){2R\ipR) (5a) 

- E ^R^R'PR + E ^R'^RPR' ■ (5b) 

R' R' 

We assume that the transitions between the configurational states of the reservoir do not depend on 
the state of the system. As neither the initial condition nor the dynamics introduces any coherence 
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between the bath macrostates, their density matrix can be written as pEnv = Hr \R) {R\ where 
^R — Tys{Pr) are the populations. The fluctuations between the configurational states of the 
reservoir are, therefore, ruled by a classical evolution: dt^R = — L/^/ ^r^r' ^r + Hr' ^r'^r^r'^ 
defining the kinetic dynamics of the environment. On the other hand, the system dynamics arise 
after tracing out all internal transitions between the configurational states of the reservoir, p = 
Y.rPR' Therefore, the evolution of the system density matrix is non-Markovian, dtP = L(p) + 
jQdTj(f{t — t)p(t), where L represents the unitary dynamics in the absence of the reservoir and 
the memory kernel introduced by it. One can compute the kernel analytically when the number 
of /?-states is small. In this case, the steady state (t = 0) correlator is G^^^ (t) = Y.r,r' (0)o>^(t)) 
where Gr = \0R) {IR\. Introducing an auxiliary matrix Gr that follows the same master equation 
as the conditional density matrix Pr but with initial conditions (i. c.) that depend on a different 
state of the reservoir 7?^ i.e., {k\6r{0)\l) = 5r^]^f5ifl{k\p]^f\l), each correlator (o'^/(0)o>^(t)) can be 
computed in terms of the elements of Gr according to ((7^/(0)o>^(t)) = Y,j^j{k\GR{T)\l) {l\GR\k) = 
(1 1 Gr{t) |0). With this, the total correlator reads: 

Its Fourier transform provides the PL spectrum 5(a), A) in the fluctuating environment. The for- 
malism is general and can be applied to arbitrary quantum optical systems where fluctuations are 
believed to play an important role.'^^'^ 

Now that we have formally included the fluctuations at a microscopic level, we need to spec- 
ify their character. They have been left completely arbitrary in the discussion until now. Since 
Gauss introduced his law of errors,^ the so-called "normal distribution" became the archetypal 
fluctuation in physical systems. Its importance and universality stems from the central limit the- 
orem, that states that the normal distribution arises from averaging a sufficiently large number of 
random variables, regardless of which distribution rules their fluctuations on an individual basis 
(provided some properties to which we return shortly). While the normal (ergodic) type of fluctu- 
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ations describes much of the physical word, it can fail dramatically for highly complex systems. 
A different class of fluctuations, for which the central limit theorem does not hold, provides a new 
paradigm. They are known as power law or fat-tail type of distributions,'^^ the most famous and 
widespread type of which is, in physics, the Lorentzian. The main property of such fluctuations is 
that they are scalefree, that is, the fluctuator can wander arbitrarily far from its most likely position, 
in contrast with a normal distribution where 5(7 deviations are so unlikely that they are regarded 
as "discovery" of a new physics (the concept of "standard deviation" is rooted in the normal dis- 
tribution and has no counterpart for scalefree fluctuators). When it is combined with impact of 
the outliers, scalefree fluctuations give rise to new notions such as kurtosis-risk, extreme-risk'^Sl 
and "black swans"'^^' — highly unlikely events that can totally alter the trajectory or response of the 




In the following, we assume that our elementary quantum optical system is embedded in such 
an environment that has scalefree fluctuations. We will assume a Lorentzian type of fluctuation, but 
qualitatively similar results would follow from other types of scalefree fluctuations]^ We are thus 
left to express the rates (l)R^Rf which enforce the steady state solution of dt^R = — (^r^r'^r + 
Hr' ^r'^r^r' to be proportional to the Lorentzian distribution ^r — ^ i {yf/2Y+{r3Y ' ^^^^ ~ 
1 . ^ is the normalization constant that removes the units (and can be simplified in the continuum 
limit) = I/? i [yp/lY^^Rdf ^ g/^^i (7^/2p+x2 = \' This is obtained with rates of the form: 

with the speed of fluctuations. They only depend on the final state. The transition towards 
the central point has the largest rate, being the most likely one no matter from which point the 
transition is initiated. 

We can now solve these equations numerically. A typical result is presented in[2f c) and, clearly, 

^As in our case the phenomenon is reproducible, it would be more accurately referred to as a "gray swan" by 
proponents of extreme outliers theories. 

^The formalism we have presented can be used for any type of fluctuations, including those of the Gaussian 
type, the case of most likely interest. Applied to the system reported in this letter, Gaussian fluctuations result in a 
broadening only of the PL line and no FIL. 
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reproduces the same phenomenology observed in our experiments; namely the unexpected emis- 
sion line, denoted "FIL" for Fluctuation Induced Luminescence. This new line arises from the 
fluctuations as sketched in[2]^d): when the system is far detuned, it infrequently wanders into reso- 
nance with the laser by the very nature of the scalefree fluctuation that allows such giant deviations. 
The photo-generated electron-hole pair will then relax into the s-shell before emitting light. De- 
pending on the speed of the fluctuations, the emitted photon can either have the frequency of the 
most likely s-shell energy when the fluctuations are faster than the recombination time, or it can 
have a frequency which is detuned from the most likely s-shell energy when the fluctuations are 
slower than the recombination time. The strong absorption gain which occurs there compensates 
for the combined scarcity and brevity of these events. This is depicted in |2]^d) with Lorentzian 
fluctuations of the levels as a function of time with three highlighted emission processes: i) from 
an extreme outlier which emits far from the most likely s-state energy, producing no notable lumi- 
nescence as such events are rare and emit weakly, ii) from a state close to the most likely s-state 
energy, producing PL as a result of the many times this weak-emission configuration is realized, 
and iii) from an extreme outlier which, by chance, brings the system into resonance, producing 
detectable luminescence from the strong emission compensating for the rarity of such events. Our 
system, therefore, implements "black-swans" in the solid state: extremely rare occurrences that 
have a huge impact and result in a qualitative change in the system. 

One can study from numerical solutions the properties of the FIL peak. Depending on param- 
eters (such as when 721 < Jf and is small), the FIL peak is sharper and brighter than the PL. 
Interestingly, it is found that the speed of fluctuations is an important ingredient that determines 
the form of the detuning dependent FIL and PL emission. Faster fluctuations decrease the intensity 
of FIL leaving its linewidth essentially unchanged. In the limit of slow fluctuations, numerical so- 
lutions converge towards the analytical expression obtained by integrating the photoluminescence 
spectrum over its fluctuations: 5(a), A) = /^oo 5'(a), 5) ^ (y^/2)H(^i-A)2 P^^Sgi^S Eq. (??) in this 
expression, one can compute the slow Lorentzian fluctuations PL spectrum in the linear regime in 
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a closed form: it is the sum of two peaks, with their Lorentzian and dispersive parts given by: 



S{(0,A) 



(3^)2+^2 (Ml)2 + ^2 

ci>FiL = , Tfil = 7io + 721 , LpiL = Yf 
(OpL = A, Tpl = no + Tf , ^PL = 721 



p=FIL,PL 



y 1 ^ Lprp/2-Kp{o)-o)p) i 
^ ;rL(7p/2)2 + (a)-a)p)2j' 

, ^nL = -7F72iA, (8b) 
^pl = 7f72iA. (8c) 



This fully describes the properties of the peaks depending on the system parameters and the ampli- 
tude of fluctuations when they are slow. This shows in particular that the peak which is sharper is 
also the one that is the more intense. Not all these attributes are conserved in the exact numerical 
solution when the speed of fluctuations is not very slow, but at least the qualitative result is spelt 
out and easily understood in this limiting case. 

The FIL peak in the experiment is observed when the excitation laser is detuned by A < 0.5meV 
from the excited orbital state resonance. At the largest detunings where FIL is still clearly observ- 
able, defined as being when its intensity is at least 1/4 of that of the PL, the system would be more 
than 20 standard deviations away from the median if it would follow a normal distribution. The 
probability for the system to fluctuate into resonance according to the Lorentzian distribution is 
approximately one in a million. However, since the quantum dynamics occurs over picosecond 
timescales, and the spectra are recorded with integration times extending over several minutes, this 
still leaves sufficient opportunities for these extremely rare events to occur and give rise to the FIL 
peak in the spectrum. They are observable as compared to off -resonant cases thanks to their huge 
impact. The same probability is completely negligible according to the Gaussian distribution and 
would not be once realized, even if running the experiment over several lifetimes of the universe. 

Fingerprints of non-ergodic fluctuations are actually commonplace in the emission from indi- 
vidual solid state emitters, including flurophores such as molecules,!^ proteins and polymers ,'^^'221 
semiconductor quantum dots,'^^'^ nanorodsl^ and nanowires. Ubiquitous phenomenology such 
as fluorescence intermittency'^^'^ and spectral wandering are observed arising from coupling to 
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Figure 3: Fluctuations due to the voltage source. The upper row shows a trace of the applied 
voltage over a 10ms window, for three sources: (a) a standard Keithley KE2400 source, (b) a 
Yukogawa source and (c) a simple battery. Beside larger amplitudes of fluctuations the Keithley 
and Yukogawa have fluctuations to all orders, that is, infrequent but large deviation from the applied 
voltage. These are the cause of FIL. 



their local environment. In our system, the source of the fluctuation was traced to be fluctuations 
of the programmable voltage source in which electronics stabilize the output voltage. Scalefree 
fluctuations are of an entirely counterintuitive character and the cause of many unpredictable and 
important events in complex environments, such as catastrophes on financial markets or sudden 
large-scale changes in meteorological systems. In our case, tracing their origin to the electronics 
in the setup was a difficult and time-consuming task since such fluctuations are typically invisible 
to the spectral noise of the devices. Furthermore, such unlikely events do not appear in the techni- 
cal description of the source provided by the manufacturer. It is not established, for instance, that 
they are exactly of the Lorentzian type, but their fat-tail property is mandatory to understand our 
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experimental observations. In most experiments, such extremely rare fluctuations of voltage have 
no impact whatsoever. It requires the black- swan character, only met in particular configurations 
such as ours, stemming from the extremely large gain of the QD absorption driven quasi resonantly 
that make these rare events impactful. In our QDs and QDMs, their effect was remarkable since 
they produced a new line in the emission spectrum. The FIL shown in [2] was obtained with a 
Keithley 2400 voltage source. We have repeated the experiment using a Yokogawa voltage source, 
and observed that the FIL peak reduced significantly in intensity, although it did not disappear 
completely. We have been able to remove the scalefree fluctuations entirely by using a simple 
battery-driven voltage follower circuit. The voltage traces measured for these three sources are 
shown in[3j The KE2400 series clearly exhibits random fluctuations with large deviations from the 
median, resulting in a clear characteristic fat-tail profile in its distribution (see the lower panel in 
[3]). The Yokogawa source, on the other hand, appears much more Gaussian-like with no observable 
large departures. The fact that a FIL is still observed shows that it, too, features extreme outliers 
with scalefree properties, but that we are unable to record them in the time windows over which 
the voltage was monitored. The difficulty to detect and characterize power-law distributions is a 
well-known problem of statistical analysis.!^ By turning these elusive outliers into black swans, 
our system therefore implements a way to probe and characterize such fluctuations. It could allow, 
for instance, to quantify the range over which the power law holds, one of the most difficult prac- 
tical problem in the field. Such fluctuations could also be intrinsic, for instance from trapping of 
a carrier in a nearby defect, with fluctuations in the induced electric field correlated with the time 
spent by the carrier in the trap, which being memoryless in time would result in Lorentzian fluc- 
tuations. In our case, however, we can rule out such possibilities from a temperature series for the 
conventional and fluctuation-induced luminescence, both shown in |4| While the conventional PL 
exhibits the characteristic phonon- sidebands with their hallmark asymmetry at low temperature,'^ 
the FIL is temperature independent, showing that it has no connection with the intrinsic fluctuation 
spectrum or dynamic of the system. 

In conclusion, we have presented a general theory of fluctuations in a dissipative quantum 
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Figure 4: Temperature series for the cross-section of (a) luminescence at resonance and (c) FIL. 
While the former exhibits a clear and characteristic temperature dependence, the FIL is temperature 
independent, pointing at an extrinsic origin for the scalefree fluctuations causing it. 

optical system and applied it to the case of a resonantly driven quantum dot in the presence of 
scalefree fluctuators. We have shown that in a system exhibiting strong gain, the combination of 
rare and impactful events may result in drastic and qualitative changes in the system, known in 
the stochastic literature as "black swan events". In our case, a new sharp line pinned at the laser 
excitation is produced. This type of luminescence has never been reported before to the best of 
our knowledge, and further investigations are needed to study whether they might prove useful for 
applications and/or characterization. The system we have presented could be used as a sensitive 
probe of its fluctuating environment and a powerful gauge of power-law tails. Although such 
fluctuations could be intrinsic to the system, we have reported here an experimental occurrence 
that was put forward, most unexpectedly, by the voltage source apparatus, and an enduring puzzle 
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was put to rest with a simple battery. 
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